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Abstract 

We discuss inelastic collisions of two rotating disks by using the conservation laws for baryonic 
mass and angular momentum. In particular, we formulate conditions for the formation of a new 
disk after the collision and calculate the total energy loss to obtain upper limits for the emitted 
gravitational energy. 
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I. INTRODUCTION 

Disk-like matter configurations play an important role in astrophysics (e.g. as models for 
galaxies, accretion disks or intermediate phases in the merger process of two neutron stars). 
The simplest models for such configurations are disks of dust. From a mathematical point 
of view they are solutions to boundary value problems of the Einstein equations. Explicit 
solutions in terms of standard functions or integrals are known for rigidly rotating disks 
of dust (the only known solutions for isolated rigidly rotating bodies) and 

counter-rotating disks of dust consisting of clockwise and counter-clockwise rotating dust 
particles 

In this paper we want to study inelastic collisions of rotating disks. (Due to the gravi- 
tational radiation there exist no elastic collisions in General Relativity.) A rigorous math- 
ematical description of such merging processes is extremely difficult, for the solution of the 
corresponding initial boundary value problem requires extensive numerical investigations. 
We adopt a different method and perform a "thermodynamic" analysis. In this way we 
forgo the detailed analysis in favour of a simpler description leading to a "rough" picture 
of the merging processes. A classical example for this procedure was given by Hawking and 
Ellis who discussed the efficiency of the collision and coalescence of two black holes, cf. Qj. 
By using the area theorem for black holes they obtained, for spherically symmetric black 
holes, an upper limit for the efficiency of conversion of mass into gravitational radiation of 
1 - I/V2 ^ 29.3%. One of our results will be a similar limit for the coalescence of two 
disks as an example for "normal matter" collisions. Considerations like these are typical 
for thermodynamics, in which initial and final equilibrium states are linked by conserved 
quantities bridging the intermediate non-equilibrium states of the system. Interestingly, a 
Gibbs equation for the thermodynamical potential energy ( -mass ) M as a function of bary- 
onic mass Mq and angular momentum J can be formulated even in the case of rotating dust 
matter, see Eq. (fTH|) . 

In particular, we study the "head on" collision of two aligned rigidly rotating disks of 
dust with parallel [scenario (a)] or antiparallel [scenario (b)] angular momenta, cf. Fig. ^ 
Rigid rotation is a universal limit for rotating disks of dust. Any amount of friction between 
the rings comprising the disk of dust will lead to an equilibrium state with constant angular 
velocity Q after a sufficiently long time. As a consequence, rigidly rotating disks are char- 
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FIG. 1: Illustration of the two collision scenarios 



acterized by an extremum in the binding energy, compared to differentially rotating disks 
with the same baryonic mass and angular momentum, see appendix El We assume that the 
initial distance between the disks is large enough to keep the initial gravitational interaction 
very small. 

As already mentioned, the dynamics of the collision process is outside the scope of our 
considerations. However, we know that the total baryonic mass Mq and the total angular 
momentum J are conserved. Finally, due to the outgoing gravitational radiation and possible 
dissipative processes, the collision ends in a stationary (and axisymmetric) configuration with 
the total baryonic mass Mq and the total angular momentum J. 

In this paper we confine ourselves to two problems: the formation of a rigidly rotating 
disk (RR disk) from two rigidly rotating initial disks and, as a second example, the formation 
of the rigidly counter-rotating disk described by Bardeen, and Morgan and Morgan 0] (RCR 
disk) from two rigidly rotating initial disks with opposite angular momenta. We will discuss 
questions like these: For which parameter values of the initial disks can the collision lead 
to a rigidly rotating or rigidly counter-rotating disk at all? Which domain of the Mq-J- 
parameter space can be reached by such processes? As we will see in section III A 21 the 
relative binding energy of the rigidly counter-rotating disks takes positive as well as 
negative values. Therefore, it is an interesting question if the formation of RCR disks with 
a negative binding energy is a possible result of collision processes. Finally, we calculate 
upper limits for the energy loss due to gravitational radiation in such collision processes. 

The mathematical analysis of these problems requires the examination of the "thermody- 
namics" of rigidly rotating disks and rigidly counter-rotating disks and the discussion of the 
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equations of state for mass and angular momentum. This is done in section [n] In section UTTl 
we discuss limits for the formation of disks after the collision. The conservation equations 
can be used to calculate the parameters of the merged disks in terms of elliptic functions. 
This analysis and the resulting plots can be found in section HVl The energy loss, i.e. the 
efficiency of the two scenarios is calculated in section |^ 

The metric coefficients of the rigidly rotating disk of dust solution are given in terms of 
ultraelliptic theta functions which reduce to elliptic functions along the axis of symmetry. 
Since the multipole moments for energy (-mass) M and angular momentum J, as the central 
quantities of our thermodynamic considerations, can be read off from the axis values of 
the metric, our analysis has to make extensive use of elliptic functions. To avoid lengthy 
calculations in the main body of the text, we relegate the analytic expressions to the appendix 
and present the results in graphical form in the main text. 

For the sake of simplicity, we restrict ourselves to identical initial disks (i.e. disks of equal 
baryonic mass and equal absolute values for the angular momenta). The generalization to 
different disks is a straightforward procedure. 



II. "THERMODYNAMICS" OF DISK MODELS 



Disk of dust solutions 



1. Rigidly rotating disks of dust 



This section is devoted to a thermodynamic description of the "ingredients" of the col- 
lision processes: rigidly rotating disks (RR disks) and rigidly counter-rotating disks (RCR 
disks) . 

The free boundary value problem for the relativistic rigidly rotating disk of dust was 
approximately discussed by Bardeen and Wagoner and analytically solved in terms of 
ultraelliptic theta functions by Neugebauer and Meinel J SQl using the Inverse Scattering 
Method. For a discussion of the physical properties see [10|. The solution is stationary 
(Killing vector: ^*) and axisymmetric (Killing vector: r/*). Its line element together with the 
Killing vectors can therefore be written in the Weyl-Lewis-Papapetrou standard form 



ds^ = e-2^[e2^(dp2 + dC') + pM^'] - e2^(dt + ady^f, C = = K 
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where U, k and a are functions of p and ( alone and is the four-dimensional Kronecker 
symbol. Note that we use the normalized units where c = 1 for the speed of light and 
G = 1 for Newton's gravitational constant. The solution can be written in terms of the 
complex Ernst potential /(p, Q = e^'-^^''''^^ + ib{p, (), where the imaginary part is related to a 
by a^p = pe~^^6^(j and a^(^ = —pe~^^b^p. In this formulation, the vacuum Einstein equations 
are equivalent to the Ernst equation 

(3f?/)(/,pp + /,cc + ^/,p) = /,P + /,?- (2) 

{k can be calculated via a path integral from the Ernst potential /.) 

The matter of the disk of dust is described by the energy-momentum tensor 

r^=£«V, £ = a(p)5(C), (3) 

where e, a and are the mass density, the surface mass density (cr(p) = if p > po? Po being 
the coordinate radius of the disk) and the four-velocity of the dust particles, respectively. For 
rigidly rotating bodies, the four-velocity is a linear combination of the two killing vectors, 

«' = e-^°(e + fir/0, (4) 

where VL is the constant angular velocity of the body. Here, as a consequence of the geodesic 
motion of the dust particles, the coefficient e"^" turns out to be a constant too, 

Vo = constant. (5) 

The RR disk solution depends on two parameters. As an example, one may choose the 
coefficients e~^° and f2e~^° of the linear combination ^ or, alternatively the coordinate 
radius po of the disk and a centrifugal parameter p = 2VL^ p'^e~'^^'^ (p ^ turns out to be the 
Newtonian limit and p — 4.62966 . . . the ultrarelativistic limit, where the disk approaches 
the extreme Kerr black hole, cf. 3 and for these and further properties). 

The baryonic mass Mq, the gravitational (ADM) mass M and the angular momentum J 
of the disk are given by 

Mo= e^f^u^d^x, (6) 



E 



M = 2 I (T,^- ^Tgij ] nX'dV, (7) 



E 
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J = - y T,,n'r]'dV, (8) 
s 

with Tij as in (jS)). S is the spacehke hypersurface t = constant with the unit future-pointing 
normal vector n*. 

2. Rigidly counter-rotating disks of dust 

An interesting example of a counter-rotating disk is the RCR disk by Bardeen, Morgan 
and Morgan j^, consisting of a clockwise and a counter-clockwise rotating component of 
dust. All mass elements move along geodesic lines and the two components have constant 
angular velocities with opposite signs. Since the net angular momentum of the disk vanishes, 
its metric can be written in the static Weyl-Lewis-Papapetrou form [a = in ((H)], 

d,2 = e-2^[e2^(d/ + da+/d^Ve'W, C = 51 r/^ = 5;, fr^. = 0. (9) 

The metric functions U and k are standard integrals, see appendix |Dl Here the energy- 
momentum tensor is the superposition of two expressions Q, 

r^' = le(uV + wV), u' = e-^^iC + Vtrf), v' = e-^«(f - Qt]'), (10) 

where and f * denote the four-velocities of the counter-clockwise and clockwise moving 
dust particles. Just as for the case of the RR disk, Vq turns out to be a constant as a 
consequence of the constant angular velocity fl and the geodesic motion of the dust particles, 
and the solution is again governed by two parameters. As with the RR disks, these could 
be chosen to be the coordinate radius po of the disk and the centrifugal parameter fi = 
2fi^PQe~^^°. However, it turns out that, instead of p, the parameter b = fipoe~^^° simplifies 
the discussions. (6 ^ is the Newtonian limit and b —>■ oo the ultrarelativistic limit.) 

The baryonic mass Mq and the gravitational mass M of the RCR disk and the angular 
momenta ±J of the counter-clockwise and clockwise rotating part (the resulting angular 
momentum vanishes) can again be calculated from the Eqs. ©-(0) (where in the formula 
for J only the energy-momentum tensor of the counter-clockwise rotating dust component 
is used). For the calculation we refer to appendix O 
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B. Equilibrium and stability 



Equili 

HQ). 



Equilibrium configurations can be described with the aid of variational principles (cf. 
For disks of dust we may consider the thermodynamic potential 

E--= J (^■^ + e^V^d'x + nnj+^, (11) 

t=to 

where R is the Ricci scalar and n indicates the number of dust components (n = 1 for RR 
disks and n = 2 for RCR disks) . The variation of E leads to 

5E = -^ j (^R'^ - ^g'^ + SttT^^ V^Sg^j d^x + e^^SMo + nnSJ, (12) 
with T*-' from (jH)) or (jl(J|) and Mq, M and J from ©-(jH))- Obviously, for fixed baryonic mass 



\Mo,J 



Mq and fixed angular momentum J, i.e. 6Mq = and 6 J = 0, the condition 6E\ 
leads to the Einsteinian field equations. On the other hand, for solutions to the field 
equations it turns out that E is equal to the gravitational mass M. Hence, one obtains 

6M = e^°6Mo + nn6J. (13) 

To illustrate the meaning of the potential E we should mention that the Newtonian limit 
of Eq. (fTTj) is given by 

E = Mo + ll p^f/^d3a: + n^, (14) 

where , U'^ and / are the Newtonian mass density, the Newtonian gravitational potential 
and the moment of inertia of the disk, respectively. In this limit, E is the sum of the rest 
energy MqC^ (c = 1), the potential energy and the rotational energy. The function E in 
Eq. (fT^ is precisely the quantity that was used by Katz in 13] to study equilibrium and 
stability of Maclaurin and Jacobi ellipsoids in Newtonian theory. In subsection III C 21 we will 
use the relativistic generalisation (jlip of E to investigate the stability of counter-rotating 
(RCR) disks of dust. 
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C. Applications 

1. Rigidly rotating disks of dust 

The extensive quantities Mq, M and J are not independent of each other. Due to Eq. (|T!^ 
(with n = 1) they have to satisfy the Gibbs formula (see also 0, [l^ or jl^ ) 

dM = e^MMo + ndJ, (15) 

where M = M{Mq, J) is a potential in Mq, J. Other potentials can be obtained via Legendre 
transformations. As an example we may use, as a consequence of (jHI)-®, the parameter 
relation 

M = e^°Mo + 2nj (16) 
to eliminate M in ()15|1 . We arrive at 

d(fiJ) = -Jdn - e^°ModVo, (17) 

and can now use the potential O J to calculate J and M as functions of the pair Q and Vq or, 
alternatively, and fi. To get explicit expressions for these "equations of state" we make 
use of the disk of dust solution |l[ Q, 0] to obtain 

fij = — y e-'^o'^^ix) dx - — 6o(/i), (18) 



where bo is the imaginary part of the Ernst potential in the center of the disk, bo = b[p = 
0, C = 0"*"), and, as a consequence of (fT7|) . 



Moin,,) = -e--^'-^ 



dVo 

d{QJ) 



J e-^"(^)6[,(x)dx, (19) 



e^» 







The parameter relation ()16p takes the form 



I e-^°(^)6[,(x)dx-^6o(/^). (20) 





M = e^oMo + 21]J=-e^oMo-^6o(/i). (21) 

For explicit expressions for the metric coefficients &o(/^)) Vo(/i), the masses Mq and M and 
the angular momentum J in terms of elliptic functions see appendix iBl 
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FIG. 2: JIMq, Mq/J and the radii po and pc as functions of the relative binding energy = 
1 — M/Mq for the rigidly rotating disk of dust. 

There is an interesting scaling behaviour of the disk parameters. For example ^Ipo, flM, 
QMq, Q^J, M/Mq, M'^/J and Mq/J depend only on the centrifugal parameter fj, but not 
on a second parameter, cf. (fTUjl - ipTj) and appendix iBl 

For an illustration of the equations of state Fig. |21 shows relations between QMq, Mq/J 
and the relative binding energy Et^ = 1 — M/Mq. The figure also displays the coordi- 
nate radius po and the "circumferential radius" pc (defined as pc = ^ J ds\p=pg^t=to,(=o = 
y/g^\p=po,t=to,c=o)- Thereby i?b — ^ and E]^ — » 0.3733. . . are the Newtonian and the ul- 
trarelativistic limit, respectively. The picture demonstrates the "parametric collapse" of 
a disk towards the black hole limit Consider a disk with a fixed number of baryons 

(i.e. fixed Mo) occupying states with decreasing energy M. Then, the angular velocity 
increases and the angular momentum decreases while the disk shrinks (po ~^ 0, while the 
"true" radius remains strictly positive, pc 1.4372 . . . ■ Mq). In the limit i?b = 0.3733 . . . , 
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2. Rigidly counter-rotating disks of dust 

As in the case of RR disks we can formulate a Gibbs relation for RCR disks, too. From 
Eq. (fT!?jl (with n = 2) we obtain 

dM = e^MMo + 2nd J. (22) 
A Legendre transformation leads to the potential QJ satisfying the equation 

d(OJ) = -Jdn - ^e^°ModVo, (23) 
From the analytic solution we find 

^ -,r^ 1 /" / 1 9v„\ tarctant , 

Vo) = — - / , - e^^° — dt. (24) 



This potential can be used to calculate the baryonic mass Mq and the angular momentum 
J via 

, M = e^«Mo + 4fiJ, (25) 



Mo ^ -2e-o ^ 



d{QJ) 

u — 



Vo 

where Vq and fl are related to the parameters po aiid b by the equations 

g-4yo = 1 + 4^2^ _ 6/V1 + 462. (26) 

The RCR disks show the same scaling behaviour of the disk parameters as the RR disks, 
e.g. Qpo, QM, QMo, Q?J, M/Mq, M^/J and M^/J depend only on the parameter b and 
not on the coordinate radius po- 

An interesting feature of the RCR disks is the strange behaviour of the binding energy, 
which is negative for b > 4.1074... as shown in Fig. El Interestingly, there are several 
physical arguments against the formation of RCR disks in this region. In fact, the following 
application of the stability analysis of Sec. Ill Bl will show a transition from stability to 
instability at 6 = 1.3393 . . . , even before the binding energy E}^ becomes negative. Moreover, 
using a general argument resulting from the baryonic mass conservation we will show in the 
next section that RCR disks with negative binding energy cannot form by collision. 
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FIG. 3: Left picture: The relative binding energy Ey, = 1 — M/Mq of the RCR disk as a function 
of the centrifugal parameter b. E\^ reaches negative values in the relativistic region (large b). Right 
picture: The coordinate radius po and the "circumferential radius" pc divided by the baryonic mass 
Mo as functions of the centrifugal parameter b. 

Fig. El also shows Pq/Mq and Pc/Mq as functions of b (right picture), with pc defined as 
in subsection III C II For increasing b the radii first decrease, reach a minimum and increase 
again. There is no "parametric collapse" of the RCR disks towards a black hole. 

In order to understand the strange behaviour of the binding energy we now apply the 
stability analysis of Sec. Ill BI to RCR disks. Calculating the Ricci scalar R of the metric Q 
we obtain from 

" 1 



E 



in 



dpdv^dC + 2nj. 



(27) 



Note that the term M/2 in (lllj) is compensated by a surface term resulting from the inte- 



gration over R. The variation 6 El 



Mo,J 



leads to the equations 



AU = 47r5(p)<5(C), S(p)5(C) := -e^'^-^^iT^ - T^) = ee'' 



2 J2-4:U\ 



(28) 



Mn 



ee 



-(7 dpdv^dC = constant, J 



■ -V i 

-e rjiU 



dpd(y9d^ = constant. (29) 



Using the relation (|28|) and the equation e^^ = e^^(l — Vt^p^e ^^) as a consequence of 
u'Ui = = — 1, (|77|) can be rewritten as 



E = 27T 



AU 



n^p^ 



■S 



dp + 2nj. 



(30) 



C Applications 



12 



In the same way we calculate the baryonic mass Mq and the angular momentum J as defined 
in (j2ni), 

PO , 



Mo = 2n I , nl. Sip)e^pdp, (31) 





PO 







Now we consider the two-parametric family of functions 

S* (p; po,b) = — , arctan , ^" (33) 

and the corresponding disk potential [as a solution of 



>p2 



1 1+ /I + 462' 
f/(p; po, &) = - In ^, \ (34) 

where 6 and po are arbitrary constants. Obviously, E depends on the parameters b, po and 
Q, E = E{b,po,Q). It should be emphasized that S{p;pQ,b) and U{p;po,b) define a family 
of trial functions which do not satisfy the equation ((201) • Only if b, po and Q are connected 
by the relation ()26|) we arrive at the RCR disk solution. To calculate the extremum of E for 
fixed values of Mq, J {SE^^^^^^ ^ = 0) we replace po and Vl via (jHTj) . (jH^ by Mq and J. With 
the explicit expressions ((HHj) and ((Mj). Mq, J and take the form 

Mo = po^i(6,^^Po), J = nplg2{b,npo), E = pog3{b,Qpo) + 2QJ (35) 

with functions gi, g2 and (73 expressible in terms of integrals resulting from ()H()j) - ()H2j) . The 
combination 

.:^Af„VJ= „'"<''-,%°' (36) 

allows to express Qpo in terms of b and s, f2po = (74(6; s). Hence, we finally obtain E in the 
form 

E = Mo^^[^'^^(^'^)] ^|f^[^'^^(^'^)] =: -MoE(6; .), (37) 

gi[b, gA{b;s)\ 

where the dependence of -E on 6 and s is given implicitly by some integral relations. The 
minima of E for fixed Mq and J are the maxima of E for fixed s and vice versa. 

The numerical discussion of E{b; s) leads to the following results: The extremum condition 
dE/db = yields the parameter relation ^E^. That means the RCR disk solution is an 
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FIG. 4: The pair of conjugate variables s and K{s) for the RCR disk. At the point A the parameter 
s takes the "critical" value s = 0.9634 . . . (corresponding to 6 = 1.3393 . . . ). The vertical tangent 
at that point indicates a change of stability. 

extremum of E (at least a stationary point). As it was shown by Poincare and discussed by 
Katz in stability can be analysed with the help of conjugate variables with respect to a 
thermodynamic potential. Here, we may choose the variable s and its conjugate K{s) with 
respect to E, 

Kis):=—[K{s),s], (38) 

where b = bf,{s) is the equilibrium relation between b and s. Fig. H] shows the pair [s, -ft'(s)]. 
The criterion for a change of stability is the existence of a vertical tangent to this curve. 
Obviously, such a behaviour is given at the point A where s = 0.9634 ... (6 = 1.3393 . . . ). 
A careful numerical analysis of the curve K{s) shows that there is no other point with a 
vertical tangent. Since a stable branch can be identified as the one with a positive slope 
near the vertical tangent, we conclude that the lower branch in Fig. |3 is stable according 
to Poincare's definition and the upper branch (dashed curve) is unstable: RCR disks are 
unstable for b > 1.3393 . . . This result includes the instabilty of RCR disks with negative 
binding energy Eb (-Eb < for b > 4.1074 . . . ). 

Because of the definition (jHUj) . Mq/J takes its maximum just at the critical value b = 
1.3393... (see Fig.^. 



A Conservation laws 



14 



III. PARTICULAR COLLISION PROCESSES 
A. Conservation laws 

As mentioned before, we consider "head on" collisions of rigidly rotating disks of dust 
for parallel and antiparallel angular momenta as sketched in Fig. ^ [from now on denoted as 
scenario (a) and scenario (b), respectively]. To compare the initial disks with the resulting 
merged disk we will make use of conservation laws. Obviously, there are two conserved 
quantities. One of them is the baryonic mass. Considering two colliding bodies A and B, 
we have 

Mo = + M^, (39) 

where the baryonic mass Mq of the final body is the sum of the baryonic masses and 
Mq of the colliding bodies (from now on, tildes denote quantities of the final bodies). A 
first consequence of (jH^ is that the (inelastic) collision of two bodies with positive binding 
energies cannot lead to a body with a negative binding energy. Namely, 

Mq-M = (M^ - M^) + (Mg^ - M^) + (M^ + -M). (40) 

According to our assumption, the first two terms on the right hand side are positive. Due 
to the loss of energy by gravitational radiation the last term has to be positive, too, 

M < + M^. (41) 

Hence {Mq > 0) the relative binding energy of the resulting body is positive, 

M 

= 1 - ^ > 0. (42) 

Mo 

Applying this result to our disks of dust we may exclude the "strange" branch (Fig. ^ 
b > 4.1074 . . . ) of the RCR disk solution: RCR disks with b > 4.1074 . . . , < cannot 
form by collisions. 

From now on, we confine ourselves to colliding RR disks with equal baryonic masses, 
= M^ = Mq. Then, Eq. §^ takes the form 

Mo = 2Mo. (43) 

(For dust, the conservation of the baryonic mass is a consequence of the local energy- 
momentum conservation T^-^-j = 0.) 
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Due to the existence of an azimuthal killing vector 77* the angular momentum is conserved 
as well. In the first scenario, the two parallel angular momenta J of the two rigidly rotating 
disks sum up to 2J in the final disk, 

J = 2 J. (44) 

In a second case we will study questions connected with the formation of the RCR disk 
by collisions of two one-component disks with antiparallel angular momenta as described in 
section III A II and assume the conservation of the angular momentum for each component 
separately, 

± J = ± J (45) 

(the resulting angular momentum of the initial RR disks and the final RCR disk vanishes). 
It should be emphasized that the assumption is not very realistic. A small amount of 
friction between the two dust components would violate the separate conservation of the 
angular momenta. Nevertheless an assumption like ()45|) can lead to deeper insight into the 
physical processes connected with the formation of counter-rotating disks and allows, by way 
of example, a comparison between collisions of disks with parallel and antiparallel angular 
momenta. 

The conservation equations enable us to calculate the parameters characterizing 

the final disk as functions of the parameters of the initial RR disks without studying the 
intermediate, complicated dynamical process. 



B. Collision of disks with parallel angular momenta 

The equations in section III A II show that the combination Mq / J of the conserved quan- 
tities depends on /z alone (and not on some second parameter). Therefore the relation 

J- ^ ^{f.) = 2^(/.) (46) 

combining the equations ()43|) and (jl^ allows us to calculate the parameter jl of the final disk 
as a function of yU alone. Using the formulae ()B2|1 and ()B3|) of appendix iBl Fig. 1^ shows that 
the initial disk ratio Mq / J as a function of yU reaches a maximum of max (Mq / J) = 2.5460 . . . 
at the ultrarelativistic limit of the disk of dust /i = 4.62966 . . . Thus Eq. ()46j) cannot have 
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FIG. 5: Mg/J of the rigidly rotating disk as function of the centrifugal parameter n 
solutions for all values of /i. Obviously, the solutions are restricted to the interval 



cf. the firt picture of Fig. [3 

Only for initial disks in this parameter range can the collision again lead to a rigidly 
rotating disk of dust. Beyond the limit fi = 1.9540 ... a collision must lead to other final 
states, e.g. black holes or black holes surrounded by matter rings. 

On the other hand, fl can take on all values in the interval [0, 4.62966 . . . ] being considered, 
i.e. there is no restriction on the parameters of the resulting RR disks. That means every 
rigidly rotating disk can be formed in a collision of two rigidly rotating initial disks. If the 
centrifugal parameter /z of the initial disks approaches the maximum /i = 1.9540 . . . , then 
the collision leads to a rigidly rotating disk with fi = 4.62966 . . . arbitrarily close to the 
extreme Kerr black hole. 

Having solved (pUj) to obtain /2 = cf. the first graph of figured one may calculate 

the new coordinate radius po from PHj) via M{fl,pQ) = 2Mo(/i, po)- The result is plotted in 
the second graph of figure 




(47) 



that is equivalent to the interval 



0<p< 1.9540 . . . , 



(48) 
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FIG. 6: Baryonic mass Mq and angular momentum J before and after the collision, cf. Eq. (|5()j) . 
C. Collision of disks with antiparallel angular momenta 

As with the rigidly rotating disk, the ratio Mq / J for the counter-rotating disk depends 
only on a centrifugal parameter, here b, and not on po- With the definition (cf. appendix 
Id) 



m - f 



In v^lp! arctan(26) - 53 f dilogi^ 



r 46 - (^2 + In ^1+^^ arctan(26) + $5 (^dilogi^) 
and with the conservation equations PHj) and one finds 



(49) 



7(&)-/(&)=49^(/^), 



(50) 



which can be used to calculate the counter-rotating disk parameter 6 as a function of the 
initial disk parameter /i. 

Similar to the case of scenario (a), by using the formulae of the appendices |B] and O 
it turns out, that (j^Uj) does not have solutions b for all values fi in the allowed range < 
fi < 4.62966 . . . Fig. El shows that f{b) reaches a maximum of /max = 0.96344 . . . for b = 
bmax = 1.33934 . . . , while the right hand side of (j^UI) grows to much larger values. Thus the 
formation of a counter-rotating disk is only possible in the small range 



M2 1 _ 
< < -/max 

i.e. from initial disks in the parameter range 



0.2408 . . . 



(51) 



< /i < 0.101003. 



(52) 
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For values of fi in this interval Eq. ()5()|1 always has two solutions: b can be smaller or larger 
than 6max, b ^ &max- Later it will turn out, that in general only the solution b < fomax is 
of physical relevance, since the formation of a counter-rotating disk system with the larger 
value of b would only be possible if external energy were put into the system. 
Finally, if b has been calculated, can be used to determine po. 

IV. PROPERTIES OF THE MERGED DISKS 
A. General parameter conditions 

In the last section we showed which initial parameters can lead to rigidly rotating or 
counter-rotating disks after a collision. Now we will discuss the parameters of the resulting 
disks as functions of the initial parameters. For this purpose we make use of the explicit 
expressions for mass M, baryonic mass Mq and angular momentum J, (jBH) - (jB3|) and (jC6|l - 
()(^8|1 . and plug them into the conservation laws (jKHjl and (jlHjl . In this way we obtain 
the desired parameter relations implicitly in terms of elliptic functions. In Fig. [3 one can 
see a number of parameter relations generated by a numerical evaluation of those implicit 
relations (solid lines). They may be compared with the corresponding post-Newtonian 
approximations derived in subsection IIV Bl (dotted lines). As expected, the curves agree 
well for small values of the centrifugal parameter /i, but differ for larger /i even qualitatively. 

As mentioned earlier [Eq. ()48j) ]. in scenario (a) the formation of an RR disk after the 
collision is only possible for < /i < 1.9540 . . . while for the final disk all parameters 
< /t < 4.62966. . . are allowed (Fig. [7|A.). The new coordinate radius po is one half of 
Po in the Newtonian limit, but smaller in the general case. In the limit p = 1.9540... 
the collision leads to a disk in transition to the extreme Kerr black hole with vanishing 
coordinate radius, po = (Fig. 03). The angular velocity increases, maximally by a factor 
of four in the Newtonian limit (Fig. EP)- For the maximum amount of lost energy M — 2M 
during the merger process due to gravitational waves one obtains a limit of 23.8% of the total 
initial mass (Fig. [7p). 

In scenario (b) the formation of an RCR disk is only possible in the small parameter 
range < p < 0.101003 ... As mentioned before (Sec. IIII C|) . there always exist two values 
of 6, b 

^ ^max — 1.33934. . ., for a given p (or Mq/ J^. In the first case (b < femax)? the 
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scenario (a) 





FIG. 7: Parameters of the resulting disks after collision (tilded quantities) as functions of the initial 
disk parameter fi in the allowed range < /i < 1.9540 . . . [scenario (a)] and < /i < 0.101003 . . . 
[scenario (b)] compared with the post-Newtonian approximations for small ^ (dashed lines), cf. 
section CYBI u, ^2, on and M belone to the initial disks. 
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FIG. 8: Behaviour of the parameters of the resulting RCR disks [scenario (b)] for the solutions 

b > branx- 

coordinate radius po shrinks to one half of the initial radius or smaller, but never reaches 
zero (Fig. EP). The angular velocity increases by a factor of four or larger (Fig. [ZP) and the 
limit for the relative energy loss is 4.2% (Fig. EH)- In the second case (6 > &max)) the new 
coordinate radius po is arbitrarily small for small p (Fig.|H(3) and the angular velocity grows 
to infinity for p ^ (Fig. |HP). 

We have seen that the emission condition (|4H) 



2M 



(53) 



forbids the formation of RCR disks with negative binding energy (6 > 4.1074. . .). The 
preceding analysis together with the exclusion of energy transfer into the system (jSHj) enlarges 
the forbidden 6-interval and restricts the formation of RCR disks aditionally. Figs. [33 
and IHP show that Eq. (j53j) is completely satisfied for the first interval (0 < 6 < fomax = 
1.33934...) and holds for the small piece 0.0728... < p < 0.101003... of the second 
intervall corresponding to 1.33934 ...<&< 3.8038 . . . but is violated for < p < 0.0728 . . . 
which corresponds to 3.8038 ... < 6 < oo. Obviously, this implies our former result that 
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RCR disks formed by collisions have always a positive binding energy, i^b > for 6 < 
4.1074. . . 

Adding now the result of section Hi Bl where we showed the instability of RCR disks in the 
interval b > 1.3393 ... we conclude that for physical reasons we cannot expect the formation 
of RCR disks in the second branch b > 1.3393 . . . 

B. Analytic parameter conditions for weakly relativistic disks 

For most astrophysical applications, rigidly rotating disk models are characterized by 
very small values of the centrifugal parameter, /i <^ 1. For such weakly relativistic disks it 
is useful and possible to obtain analytic parameter conditions derived by post-Newtonian 
expansions of the equations of state and the conservation laws fj46|) and ()50p. (For counter- 
rotating disks, for which Eq. ()50|) has two solution branches, we restrict ourselves to solutions 
with b ^ 6inax excluding the "strange" branch with negative binding energy.) 

1. Scenario (a) 

The expansion of the relations in appendix ^ in a power series for n leads to 

^M=^ (^/^^/^ - ^/x^/^) + 0{^,y% (54) 

fiMo = ^(/.^/^-i/.^/^)+0(/.^/^), (55) 

and 

fiV=^(/.^/^-i//^)+0(/.^/^). (56) 
The conservation equation ()4(ij) now can be written as 

lO/il/2 ^ -3/2 ^ -5/2) ^ 2(10/^1/2 ^ ^3/2) ^ (^(^5/2) (57) 

with the solution 

12 

/i = 4/x--/i2 + 0(^3)_ (5g) 

With yu and the expansions and u = Qpo = 1/^2 (/i^/^ — l/2yU^/^) + 0(/i^/^) one can 
calculate po from Eq. (j^ . 

Po/po = ^-^/x + 0(p2). (59) 
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That means, in the lowest order, the new coordinate radius is one half of the original radius. 
For the change of the angular velocity, we obtain 

n/n = A-6i2 + 0{iJ^), (60) 

i.e. the rotation of the final disk will be four times faster than the rotation of the initial 
disks in the lowest order. By comparing the gravitational masses of the two initial disks and 
the final disk, one can determine the energy loss due to gravitational radiation during the 
dynamical process. The relative energy loss is 

{M -2M)/2M = -^fi + 0{fi^), (61) 

where the minus sign shows, that the energy indeed leaves the system. 

To compare the post-Newtonian approximations with the exact parameter conditions of 
the last section they are plotted in figure [7| as dotted curves. 



2. Scenario (b) 

With the expansion of Eq. pUj) in terms of the centrifugal parameter b, f{h) = ^(56 — 
fb^) + 0{P), Eq. dnOl) takes the form 

(lO/i^/^ ^ ^3/2) ^ 0(^5/2) ^ ^(5^ _ ^^3) ^ (>^(^5) (g2) 



6=v^(//^ + ^/.^/^)+0(/.^/^). (63) 



with the solution 



From this and with Mq = Apo/Sn [b^ - fb^ + 0{b^)j one finds for the coordinate radius of 
the disk 

1 197 

Po/Po = ^-Y^fi + Oifi'). (64) 
The change of the angular velocity can be calculated from Q = uj/po and fl = 



6/(poV 1 + 462). igads to 

n/n = 4 + ^fi + 0{fi^). (65) 
For the relative energy loss one obtains 

(M-2M)/2M = -^/i + C»(/i2). (66) 
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Comparing the collisions of disks with parallel and antiparallel angular momenta, it 
turns out, that to the lowest order, the change of the angular velocity and the change of 
the coordinate radius are equal in both scenarios. Namely, the lowest order terms represent 
the Newtonian results and in the Newtonian theory there is no influence of the direction of 
rotation on the gravitational field. Interestingly, a similar effect occurs in the post-Newtonian 
regime: The energy loss for both scenarios also coincides to the lowest order. 

As an example, one can calculate the energy loss for the merger of two disks with mass 
and radius of our milky way (po = 15000 pc and M = 2 • IO^^Mq, respectively). The 
corresponding centrifugal parameter /i = 3 ■ 10^^ is so small that the post-Newtonian ap- 
proximation is sufficient to calculate the efficiency. With Eq. fl(j(jj) one finds for the energy 
loss AM = M -2M = -A - lO^M^ for both scenarios. 



V. THE EFFICIENCY OF GRAVITATIONAL EMISSION 

An important quantity in the collision process is its efficiency rj, which is the negative 
value of the relative energy loss, 

' 2M ^ ^ 

Obviously, t] is the upper limit for the energy transported away by gravitational radiation. 
(Some part of the energy will also dissipate due to friction. To guarantee the formation of a 
new rigidly rotating or rigidly counter-rotating disk, friction is necessary to force a constant 
angular velocity.) 

As shown in section IIVB[ for weakly relativistic disks the efficiency is given as 

r] O.Sfi, (68) 

cf. dnn). 

In section HV Al we already presented the relative energy loss as a function of the centrifu- 
gal parameter /i of the initial disks, cf. Fig. [7p. Here we derive a more explicit expression 
for the efficiency making use of the scaling behaviour of the disks: M/Mq and Mq/J are 
functions of the centrifugal parameter alone {fi for RR and b for RCR). Therefore M/Mq is 
a function of Mq / J, 

- = FW. (69) 
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FIG. 9: The efficiency r/ for scenario (a) (long curve) and scenario (b) (short curve) as function of 
the physical quantity Mg / J of the initial RR disks up to the natural end points of the intervals 
for which disk formation is possible. 



Using the conservation laws ()43|) and ()44|) or we obtain the following expressions for 
the efficiency r/ as a function of the quantity Mq / J of the initial disks, 

F{2x) 



ri{x) 



1 - 



F{x) 



for scenario (a), 



(70) 



ri{x) = 1 Sty scenario (b), 

r \x) 



(71) 



where F and F can be calculated in terms of elliptic functions from the equations of state 
of the RR disks and the RCR disks, respectively. Accordingly, rj is determined completely 
by the relative binding energy Et^ = 1 — M/Mq = 1 — F{x) of the RR or RCR disks. 

The efficiencies (fTOj) and (f7T|) are plotted in Fig. El Note that the efficiency measures the 
energy loss. The picture shows that the formation of a counter-rotating RCR disk [scenario 
(b)] is "more efficient" than the formation of a rigidly rotating disk [scenario (a)] for initial 
disks with the same amount of angular momentum and baryonic mass, but is possible only 
in a much smaller parameter range. It ends with a comparably small maximum value of 
rj ^ 4.2%. The limit rj ^ 23.8% for the collision of disks with parallel angular momenta has 
the order of magnitude of Hawking's and Ellis' limit rj ^ 29.3% for the coalescence of two 
spherically symmetric black holes. 
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VI. CONCLUSION 

In this paper we have discussed the coUision and merger processes of rigidly rotating 
(RR) disks leading again to a RR disk or a rigidly counter-rotating (RCR) disk. The 
conservation equations for mass and angular momentum showed, that these processes are 
restricted to limited parameter intervals for the initial disks: Only for rigidly rotating disks 
whose centrifugal parameter does not exceed a given maximum value do these conditions 
allow the formation of a new rigidly rotating disk, and in an even much smaller parameter 
range, the process can lead to a counter-rotating disk. (But of course, even if the balance 
equations can be satisfied in these parameter ranges, it is not clear if these processes are 
dynamically possible.) We were able to calculate the physically relevant relations between 
the parameters ("equations of state") by a numerical evaluation of the exact but implicit 
conditions resulting from the conservation equations. Explicit analytical expressions could 
be derived for the post-Newtonian domain. 

It turned out, that every RR disk can be the result of such a collision process (every point 
of the two-dimensional parameter range can be reached) , while the formation of RCR disks 
is restricted. RCR disks with negative binding energy (as described by a branch of the RCR 
solution) cannot be formed in a collision process. An interesting result is the calculated 
upper limit for the efficiency of 77 ~ 23.8% for the formation of RR disks and of 77 ~ 4.2% 
for the formation of RCR disks. 

Counter-rotating disks are interesting initial configurations for axisymmetric collapse sce- 
narios with gravitational emission. In these cases, the mathematical analysis consists in the 
discussion of initial-boundary problems for the vacuum Einstein equations. A typical exam- 
ple will be published elsewhere. 
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APPENDIX A: WHY RIGID ROTATION? 

If there is a small amount of friction between the dust rings (which can be anticipated 
for every realistic system), differentially rotating disks will develop towards an equilibrium 
state with rigid rotation. Due to the loss of energy via gravitational or thermal radiation 
the gravitational mass M of the disk decreases in these processes. Therefore we expect to 
find rigidly rotating disks as minima of the gravitational mass M compared to differentially 
rotating disks with the same baryonic mass Mq and the same angular momentum J. Here, 
we show that M indeed takes an extremum. ^ 

Following the calculations in , one can generalize the Gibbs formula ()15j) to disks with 
differential rotation Vt = f2(p). The result is 

SM= j [n^5{pM,V^) + ^5{pj^)]d?x, (Al) 

t=to 

where /iq = /io(p) is the chemical potential (specific free enthalpy), g is the determinant 
of the metric tensor and puo and pj denote the densities of baryonic mass and angular 
momentum, respectively, defined by 



Mo= J pMoV^d'x, J = J pj^/^d-'x. (A2) 

t=to t=to 

If one now considers the set of all disks with different rotation laws Q = fi(p), but fixed 
baryonic mass and angular momentum, i.e. with the constraints 

Ci := Mo = constant, C2 '■= J = constant, (A3) 

then one can find the disks with extremal gravitational mass M (and accordingly — due to 
Mq = constant — with extremal binding energy Mq — M) . With the Lagrange multipliers 
Al and A2, the condition for a stationary point is 

6{M + XiCi + A2C2) = 0. (A4) 



^ To avoid confusion we want to point out that we use two different kinds of variational principles in this 
paper. In this section, we compare different solutions to the Einstein equations (with different rotational 
laws) to find the rigidly rotating disks as extrema of the gravitational mass M . In contrast to this we 
consider trial functions that do not solve the field equations in Sec. Ill Bl to find the solutions as extrema 
of the thermodynamic potential E. 
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Using ()A1|) - ()A3|1 . the latter equation takes the form 

J lifio + XMpmoV^) + + X2)5ipJ^/^)] d^x = 0. (A5) 

t=to 

Since pMo and pj are varied independently, one finds 

Q = constant, po = constant, (A6) 

i.e. the family of rigidly rotating disks of dust. 

The same considerations are valid for counter- rotating disks of dust. The generalization 
of Eq. (HH) is 




where the index i distinguishes between the two fluid components. With the restric- 
tion to disks with the same chemical potential {pq=Pq=: ^po) and baryonic mass density 

(1) (2) (1) (2) 

{PMo=PMo=- ^PMo) and opposite values of angular velocity {Q= — Q=: Q) and angular 
momentum density (pj= — pj='- pj), Eq. ()A7|) reduces to Eq. ()A1|) . i.e. the extremal 
problem leads to the rigidly counter-rotating disks with Vt = constant. 



APPENDIX B: MULTIPOLE MOMENTS OF THE RIGIDLY ROTATING DISK 

By specifying the formulae of the general solution to the axis of symmetry, it is possible 
to calculate all multipole moments of the disk [l^ . In particular, one can derive expressions 
for the gravitational (ADM) mass M, the baryonic mass Mq and the angular momentum J, 



nM{p) = --{uj{p)ai{p) + bo{p)) 



(Bl) 



QMoip) 



ai{p), 



(B2) 



n\j{p) 



f w(/i)ai(/i) + ho{p) 



(B3) 



where one needs the relations 



ai 



1 



2^/l + fiHoip)(h' 



m 

K{h) 



) + hip) + (1 + /x')^/^-^Ao(am(/(/i), h'), h)) 



hh' 



- iE[am{I{p), h'), h'] -h'l]+ h{p) 
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bo 



sn(/(/i), /i'(/i))dn(/(/i), h'ifi)), 



l-^cn(/(^),/.'(/i)), 



^2^o-^cn2(/(/.),/^'(/.)), 



h{fi) 



\ 



^ 1 + 



1 /■ ln(Vl + ^2 + a;) 



TT 



X 



h' = 



Vl + a/a* - ^ 



with the complete elhptic integrals = F{TT/2,k) and -^(/i;) = -E'(7r/2, fc), Heumann's 

Lambda function Ao(^/', A;) = f [£;(A;)F(7/;, A;') + K{k)E{ilj, k') - K{k)F{ilj, k')], k' = 
and the Jacobian elliptic functions sn, cn, dn and am. 



APPENDIX C: MULTIPOLE MOMENTS OF THE RIGIDLY COUNTER- 
ROTATING DISK 



For the metric by using the field equations, the relations (0)-® for mass and angular 
momentum can be expressed using only the metric potential U and its p- and (^-derivative 
in the disk. For general counter-rotating disks one finds 



PO 

M = j f/,^pdp, 





PO 



Mo = y e-^^(l - 2pf/,,)(l - pt/,p)f/,cpdp 





and 



1 



PO 



J=2l e-^%/P^,p(l-pf^,p)f^,CP'dp. 





(CI) 



(C2) 



(C3) 



These equations can be applied to rigidly counter-rotating disks (RCR disks) with the ex- 
phcit expressions (see jsj) 



f/(P,C = 0) = -ln 



2h 



1 + y/f + 462pVpg) 



(C4) 
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2^ 1 Cretan ^^£^£M 



where = 6/ (poVl + 46^) is the constant angular velocity. They lead to 



(C5) 



M = — ( 1 - — arctan(26) 



(C6) 



Po 
47r6 



1 + 462)1/4 



ln^^±i?arctan(26)-^rdnogi±^ 



and 



J 



Pi 



167r62 



Vl + 462 



46 - I 2 + In ^^^^^^ ^ arctan(26) + $5 fdilog^ ^ 



with the dilogarithm function 



(C7) 



(C8) 



dilogf^;) = / ^'^^ dw. 

1 — w 



(C9) 



APPENDIX D: THE COUNTER-ROTATING DISK OF DUST METRIC 

The Inverse Scattering Method (cf. can be used to calculate the function U{p,() 

in the line element Q of the counter-rotating RCR disk. For an angular velocity fl and a 
coordinate radius po it leads to 

1 

]-n\^ _ r^n - P^l _ _ 

(Dl) 



rr. .^ 1 f ln[l - a(l - P)] 



By using elliptical coordinates {C,,ri) with 



P = Po 



v/(i+e)(i-^2), c = Poe^, 



fD2) 



it is possible to express U in terms of the dilogarithm function ()C9|) . The result is 



U = — arctan - ■ In 

2 



L4 



1a 1 1 1 

-Artanhr? m 

2 ' 2 



~ ^ ^ + dilog AT] , 



(D3) 
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where 





1 — a 



a 



(D4) 



and 



n — 



(D5) 



The function k can be calculated from U via line integration by using the field equations 
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